General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



r 



EFFECTS OF PROCESSING ON THE CARRIER 
LIFETIME IN SILICON SOLAR CELLS 


by P. A. lies and S. I. Soclof 


GLOBE - UNION INC. 

Centralab Semiconductor 
4501 N. Arden Drive 
El Monte, California 91734 


Prepared for 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
NASA Lewis Research Center 


Contract NAS 3 - 15689 


(NASA-CR- 134523) EFFECTS OF PROCESSING OK N75-3157 

THE CARRIER LIFETIME IN SILICON SOLAR CELLS 
(Globe-Union, Inc.) 107 p EG $5.25 CSC! 10R 

Unclas 
G3/44 35276 





M 


T. Report 'No'. 

NASA CR 134520 


accession No . 


4. Title and Subtitle 

EFFECTS OF PROCESSING ON THE CARRIER 
LIFETIME IN SILICON SOLAR CELLS. 



eciprent’s catalog No. 


5 . Report Date 

February 1975 


6. Perf. Organisation Code 


7. Author (s) . 

P. A. lies and S. I. Soclof 


8. Perf. Organisation Report# 


9. performing Organisation. 

‘ Globe - Union Inc. 

Centralab Semiconductor 
4501 N. Arden Drive 
El Monte, California 91734 


10. Work Unit No. 


Contract or Grant No. 
NAS 3 - 15689 


13. Type of Report. 
Contractor Report 


14. Sponsoring Agency Code. 


12. Sponsoring Agency. 

National Aeronautics and Space Admin. 

Lewis Research Center 

21000 Brookpark Rd. Cleveland, Ohio 

44135 


5. Supplementary Notes. 

Project Manager, Henry W. Brandhorst, Jr. MS 302—1 
NASA Lewis Research Center. 

Cleveland,, Ohio 


it. Abstract. 

•Phis was a study of the effects of solar cell processing on the 
minority carrier lifetime in silicon. Measurments were made of the 
diffusion length and solar cell properties for ingots covering a 
wide range of doping and other ingot properties. The results showed 
that the cell processing did not cause appreciable change in 
diffusion length. The effects of heavy doping concentrations on 
solar cell performance were studied. 




17. Key Words (Suggested by Author (s) ) | 18 . Distribution Statement, 

j * Lifetime Unclassified - unlimited 

Silicon Solar Cells 
Process effects 


L9, Security Classif. 
Unclassified 




JO. Security Classif. (pageEl. 
Unclassified I 


ages 


* For sale by the Nat. Tech. Inf. Service, Springfield, Va. 22151 




TABLE OF CONTENTS 


Page No 


1.0 INTRODUCTION 

1.1 Objectives of the Present Study 

1.2 Dependence of Silicon Solar Cell Characteristics 
on Silicon Properties 

1.3 Current Experimental Situation 


1 

1 

2 

7 


2 i 0 EXPERIMENTAL DETAILS 

2.1 Description of Special Ingots 

2.2 Measurement Methods Used 


9 

9 

12 


3.0 EXPERIMENTAL RESULTS 

3.1 Diffusion Length Results 

3.2 Solar Cell Results 

3.3 Other Results 


16 

16 

18 

20 


4 ■, 0 DISCUSSION OF RESULTS 

4.1 Diffusion Length Measurements 

4.2 Diffusion Length Values 

4.3 Low Resistivity Silicon Results 

4.4 Relevance of Program to Objectives 

4.5 implications for Future improvements in Cell 
Output 


25 

25 

25 

26 
30 

30 


5 . 0 CONCLUSIONS 


6 . 0 RECOMMENDATIONS 


7.0 REFERENCES 


31 

32 

33 


APPENDIX A 

Discussion of Methods Considered for Measuring 
Diffusion Length 


35 


(i ) 


TABLE OF CONTENTS (contd.) 


Page No 

• -j 

APPENDIX B 

Detailed Discussion of Photovoltaic Diffusion 
Length Measurements 40 

TABLES 58 

FIGURES 70 




LIST OF TABLES 


No. Title Page No. 


1. Ingot Details 

2. Cell Processing Sequence 

3. Etch Pit Densities 

4. Saturation Current Densities and V Qc Calculations 

5. Electronvoltaic Diffusion Length Calibration 
Samples for 1^ Method 

6. Electronvoltaic Diffusion Length Calibration 
Samples for 1^ Method 

7. Measured Absorption Coefficient for Four Light 
Sources Used for PV Method 

8. Calculation of the Photon Current of an LED 

Using a Calibrated Photocell 65 

9. I, and I f L-measurements on Several Samples to 

Snow Repeat ibility 66 


10. 

I Measurements for Thick Samples from Ingot 
9/3-9 (Long Diffusion Length) 

67 

11. 

Calculated Diffusion Lengths in 
from Ingot 973-9. 

Thick 

Samples 

68 

12. 

Calculated Diffusion Lengths in 
from Ingot 973-9 

Thick 

Samples 

69 


(iii) 


58 

59 

60 
61 

62 

63 

64 


LIST OP FIGURES 


No., Title Page No. 

1. Diffusion Length vs. Resistivity (N-Si) 70 

2. Short Circuit Current Density vs. Diffusion Length 71 

3. Short Circuit Current Density vs. Resistivity 72 

4. Open Circuit Voltage vs. Resistivity (Theory and 

Experiment 73 

5. Sample Structures for PV L-measurements 74 

6. Diffusion Length after Processing Steps 

Ingots 973- (1 to 3) 75 

7. Diffusion Length after Processing Steps 

Ingots 973- (4 to 7) and 973-15 76 

8. Diffusion Length after Processing Steps 

Ingots 973- (10 to 12) 77 

9. Diffusion Length after Processing Steps 

Ingots 973- (13, 14, 16, 17) 78 

10* Diffusion Length after Various Processing Steps 79 

11. Diffusion Length vs. Resistivity 80 

12. Diffusion Length vs. Acceptor Concentration 81 

13. Minority Carrier Lifetime vs. Acceptor 

Concentration 82 

14* short Circuit Current vs. Acceptor Concentration 83 

15* Short Circuit Current vs. Diffusion Length 84 

16. Spectral Response vs. Wavelength 85 

17. Spectral Response vs. Wavelength 86 

18. Open Circuit Voltage vs. Acceptor Concentration 87 

19. Forward Voltage Distribution on Mesas for 

Three Resistivities 88 

20. Envelopes for in ! f vs. V f for Three Resistivities 89 

21. Envelopes for In I f vs. V- for Three Junction 

Depths 90 

22. Envelopes of In I- vs. V f for Three Junction 

Depths 91 

23 . Open Circuit Voltage vs. Acceptor Concentration 92 

■ dv) 


LIST OF FIGURES (contd.) 

Title 

Geometric Arrangement for Theoretical Analysis 
Later Measurement Circuit 

Front Current (Tungsten + 1.06 urn filter) vs. 
EV Diffusion Length 

Absolute Spectral Response for Cell 2A-4B 
Step Sample 

log (Back Current) vs. log (Slice Thickness) 
log (Back Current) vs. Slice Thickness 


Page No 

93 

94 

95 

96 

97 

98 

99 


1 


. , - . v - r i-T.r^.^.-.r„ ,„ -■■■* 


SVMMASY, 

This report describes measurements to determine the effects 
that solar cell processing steps have on the minority carrier 
lifetime in the silicon. 

The experimental work included measurement of diffusion 
lengths of the silicon at various stages of cell fabrication, 
of solar cell properties, and of the PN junction properties. 
The measurements covered a wide range of resistivities 
(10 to 0.01 ohm-cm) and included combinations of extreme 
concentrations (high and low) of dislocations and oxygen. 

The results showed that the major cell process steps did 
not usually cause appreciable change in diffusion lengths. 
The solar cell currents obtained supported the diffusion 
length measurements. Attempts were made to relate the 
voltage behavior of the cells with the diode properties and 
the properties of the ingots. 


EFFECTS OF PROCESSING ON THE CARRIER 
LIFETIME IN SILICON SOLAR CELLS 

1.0 INTRODUCTION 

Silicon solar cells have received intensive development, 
and they have provided secondary power for most of the 
spacecraft flown to date. The requirements for space 
useage are demanding, requiring a combination of highest 
electrical output, light weight, and a good degree of 
resistance to damage caused by charged particles encoun- 
tered in orbit. The silicon material providing the best 
combination of properties has been high quality, single 
crystal, p-type, around 1 to 10 ohm-cm, with a shallow 
N+ diffused layer. 

There is still the possibility of obtaining increased 

output from silicon solar cells. A most promising area 

would be to obtain the higher voltages theoretically 

possible when lower resistivity silicon is used as the 

(1-3) 

starting material for solar cells. 

in recent years there has been renewed interest in using 
solar cells to provide some earth surface power at levels 
sufficient to help the anticipated energy crisis, it 
is probable that in order to maintain costs sufficiently 
low, the silicon used will not be in single crystal form, 
and also may not fall in the resistivity range used 
most for space cells . 

1.1 objectives of the present Study 

This program studied solar cells made from silicon ingots 
deliberately doped to give a wide range of resistivity, 
particularly lower resistivities (below 1 ohm-cm) . 


i ..... 

t . - 

The work used present cell processing methods and evaluated 
the effect of these methods on the minority carrier 
diffusion length (or lifetime) in the silicon ingots 
and on the properties of the solar cells made by applying 
all the process steps to this range of ingots. 



(a) TO predetermined specifications, a number of special 
silicon ingots was obtained; these ingots spanned/ 
three decades of resistivity ( 0.01 to 10 ohm-cm) . 
The group also included Ingots which had the extreme 
values of two other important silicon crystal growth 
parameters, namely oxygen concentration and disloca- 
tion density. 
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(b) These special ingots were processed into solar cells . 
in addition to evaluation of cells made from this 
wide range of starting ingot properties, the minority 
carrier diffusion length was measured after various 
stages of cell processing to check if the processing 
had degraded the diffusion length of the original 
ingot. 


Dependence of Silxcon Solar ceil cnaracteriscics on 
Silicon properties 

The solar cell is a large area, shallow PN junction, 

(PN diode) formed by thermally diffusing one conductivity 
type of impurity into silicon doped with the other con- 
ductivity type impurity. For this study, the N-type 
diffused layer was formed by diffusing phosphorous into 
P-type, boron doped wafers * 

in the power generating (photovoltaic) mode of operation, 
the junction collects current carriers generated by 
absorption of sunlight at various depths of penetration 



of sunlight into the silicon. These carriers can move 
from their point of generation, and if their lifetime 
is sufficient, they move to the PN junction region and 
are collected as useful current 0 The PN junction has 
a Luilt-in voltage between the P and N regions; the 
power generated by a solar cell is derived from the 
passage of the collected charge carriers through this 
junction voltage. 

Several properties of a silicon crystal can control 
cell characteristics. 


(a) influence of Resistivity 

The resistivity of the starting silicon (P-silicon 
in this case) is determined by the boron doping 
concentration and is given by the expression 





where q is the electronic charge, 

N_ is the acceptor (boron) concentration, and 
jj, is the mobility of holes. 

The light-generated properties of the solar cell 
are related to the silicon diode characteristics 
by the following equation: 

*l - r sc - * 0 (cxp - 11 (2 > 

in this equation, the voltage V and the current I 

L ii 

represent the electrical output of an illuminated 

solar cell operating into a fixed external load* and 

I is the current obtained when the cell load termi- 
sc 

nals are short-circuited. The second term on the 
right hand side of (2) describes the pn junction 
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characteristics. Here q, k, and T have their usual 

meaning (electronic charge, Boltzmann's constant, 

and absolute temperature) , and A is a correction 

factor indicating deviation from an "ideal" diode 

(where A = 1) . The diode saturation current, i , 

o 

expresses the silicon properties which control the 

PN junction characteristics. For a simple case 

(diffusion diode theory ), i is given by: 

2 2 
n. D Vf n. D W 

*° = q n l coth r + q ri r 1 coth f- 

Dp p An n 


( 3 ) 


where ^ is the dopant concentration influencing 
* resistivity, 

L p ^ n is the minority carrier diffusion length, 
and n^, D, W refer to the intrinsic carrier concen- 
tration, the diffusion coefficient of the 
carriers, and the thickness of the silicon. 


For most solar cell structures, the base region 
component of I (the second term in (3) in this case) 
is usually very much larger than that contributed 
by the diffused layer (the first term in (3)), and 
thus the diode characteristics depend primarily on 
the silicon material properties of the p-type base; i.e., 

n. 2 D W 

_ i n . . p 

I o* q n"l coth h (3a > 

An n 

Similarly the maximum voltage obtainable from the 
solar cell at a given illumination level (corresponding 
to the case where the cell load terminals are left 
open circuited) is given by: 

V oc* IT 1W 1*°“ (4) 

o 
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High resistivities (n low) generally have high 

A 

L -values. Since I increases for higher L-values, 

SC 

high resistivity silicon can lead to higher I 

sc 

values (and conversely, low resistivity leads to 

lower I ) . 
sc' 

The obtainable range of N values (from 10 13 to 10 20 
-3 A 

cm ) is much larger than the corresponding range 

of L values (# 300 to 5 /(.m) , meaning that N has 
n A 

the strongest influence on determining the i Q -value 
in equation (3a) . 


Since high N_ means lower I , lower resistivities 
A o 

would be expected to lead to higher voltages V , 
by (4). 


The above considerations explain the statement made 

in the introduction above. This program examined 

the silicon properties and the behavior of cells 

over a fairly wide range of N and l values, to 

A n 

see if the above theoretical expressions were reli- 
able in predicting cell behavior. 


(b) Considerations on Minority Carrier Diffusion Length 
( i ) Unwanted impurities 

in addition to dependence on the background doping, 

L can be determined by the other properties of the 
silicon. 

in particular, the presence of very small concentra- 
12 14 ~3 

tions (#10 -10 x cm ) of any of several harmful 

impurities (including gold, iron, and copper) can 
cause serious decrease in L, because these impurities 
can seriously increase the carrier recombination. 
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Thus some of the precautions described below are 
empirically selected (by continued useage in normal 
cell processing) to reduce the chance of accidental 
introduction of sufficient concentration of these 
"killer" impurities. \ 

(ii) Other impurities 

Although not as serious in their effect as the impuri- 
ties discussed in (i) , there is the chance of intro- 
duction of impurities other than boron into the sili- 
con crystal. These impurities may be a byproduct 
of the purification process for the silicon used for 
the growth of the silicon crystals (e.g. carbon), 
or may be the result of the growth conditions (e.g. 
introduction of oxygen from the quartz crucibles 
used to hold the molten silicon during crystal growth) . 

(iii) Crystallographic perfection 

The diffusion length for carriers is reduced if they 
cannot move freely through the crystal lattice. 

Although all the crystals used here are single crys- 
tal, (meaning that no grain boundaries are included), 
there is still the chance of internal crystal imper- 
fections such as dislocations or clusters of impurities. 

It can now be appreciated better why the ingots specified 
for this program included a wide range of intentional 
boron doping concentrations, high and low values of oxy- 
gen concentrations (oxygen could possibly interact with 
other impurities or crystal defects, thereby affecting 
L-values), and a range of dislocation density values. 


1 


1.3 



Current Experimental Situation 

The extent of the information on the dependence of solar 
cell properties on the silicon resistivity is summarized 
by Figures 1 through 4. 

Figure 1 shows the dependence of L on resistivity for 
N-silicon; the decrease in L for decreased resistivity 
is seen. 

(4) 


Figure 2 gives the dependence of j for a cell on L; 

s c 


the measurements were made under simulated space sunlight 
(called air mass zero, AMO) . 


All cells were «*12 mils thick and were made from 10 ohm-cm 
P-silicon, and the diffusion length was systematically 
decreased from the starting values ( * 150 ^um) by exposure 
to increasing dosages of 1 Mev electrons. J is seen 

“ SC 

to decrease rapidly for L £ 100 jum; the slow rate of 
change for the higher L-values can be explained by the 
fall-off in the fraction of the solar spectrum which 
penetrates to depths greater than 150 yum; thus even if 
the carriers have greater diffusion lengths, there are 
relatively fewer carriers generated to contribute to 


I by moving to the PN junction. Using a combination 
sc 


of curves like Figures 1 and 2, the dependence of j 


sc 


on resistivity can be estimated. Figure 3 shows such 
a curve, with some measured values indicated by squares. 


Figure 4 shows the theoretical* and observed variation 
of the maximum voltage obtained from cells (under AMO 
illumination) , as the resistivity of the bulk silicon 
was varied. This figure shows again the interest in 


determining the causes for V qc fall-off for resistivities 


* using the simple theory outlined in Section 1.2 
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below 0.5 ohm-cm, because attaining V values closer 

oc 

to the theoretical values would result in significant 
increase in cell output. 


2.0 




EXPERIMENTAL DETAILS 

2.1 Description of Special ingots 

Table 1 lists the properties and specifications of the 
ingots used for this study. The ingots were boron-doped 
to be near 0.01, 0.05, 0.1, 1, and 10 ohm-cra. The ingots 
had either high or low concentrations of oxygen or dislo- 
cations. 

The oxygen content was a result of the method of ingot 
growth; ingots with high oxygen concentration were 
crucible-grown, those with low oxygen content were grown 
without a crucible. 

The dislocation density in the silicon ingots was influ- 
enced by the perfection of the seed crystal, and by the 
crystal growth conditions — mainly the speed of growth 
and the temperature gradients present while the crystal 
was growing. 

2.1.1 Sample preparation 

The ingots were all grown in the (111) direction, and 
all slices were cut perpendicular to the growth axis, 
thus ensuring that the slices were in the (111) plane. 

The slice cutting to 2 x 2 cm squares was carried out 
by an internal diameter diamond-edged saw blade, with 
water cooling to minimize crystal damage during cutting. 

For each ingot, the following slices were cut: • 

(a) A thick slice («y 0.100 in.) was cut near the end 
of each ingot and chemically polished by etching 
in a mixture of HF and HN0 3 ( ^1 HF: 4 HN0 3 ) . For 
all ingots except # 10, 11, 12, Schottky metal- 
silicon collecting barriers* were applied to one 
* see Appendix b2 for details of these barriers. 
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polished face; the opposite face was then sand- 
blasted, and an open-structure ohmic contact pat- 
tern was applied by evaporating titanium-silver 
metals through a metal shadow mask. 

Figure 5 shows the structures used. 

These samples were used to determine the diffusion 
length for the unprocessed ingot. 

(b) The remaining ingot was sliced to give slices of 
thickness greater than twice the diffusion length 
as estimated from the manufacturer's measurements 
on the ingots (or from prior experience with sili- 
con ingots) . 

These slices were all lapped to 5 jul m finish and 
were then mechanically-chemically polished, using 
the Lustrox process. 

The slices were divided into several groups. 

(i) A few slices (#5) were treated as described 
in (a) above, to provide additional samples to evalu- 
ate diffusion length before any further processing. 

2.1,2 Cell processing Sequence 

The remaining slices were taken through the cell 
processing sequence shown in .Table 2. A drop- 
sequence was used, i.e. several slices were removed 
after the major process steps shown in Table 2. 

(ii) Those slices which passed through the complete 
sequence of process steps were evaluated as solar 
cells. 
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(iii) The other slices were sandblasted on the 

back face, had an open metal contact structure applied 
to this side, and had a full contact applied to the 
diffused side (like Figure 5), for use in diffusion 
length measurements , 

(iv) Some of the slices from steps I and II were 
chemically polish-etched and were then placed in 
Sirtl etch^ 7 ^ to form etch pits at edge disloca- 
tions. For these slices the dislocation density 
was measured, to compare to the density quoted by 
the ingot supplier, and also to check whether the 
dislocation density had changed during diffusion. 

(v) On some of the diffused slices, small areas 
were isolated by masking and were then etched to 
form a number of "mesa" PN diodes. The unillumi- 
nated forward diode characteristics were measured 
on these mesa diodes. 

2.1.3 precautions Taken in the processing 

The cell process sequence used was a proven manufacturing 
procedure. As far as possible, at all stages of sample 
and cell preparation, care was taken to minimize the 
introduction of impurities which could affect the sili- 
con properties (see 1.2 (b) (i) ) . 

This care extended to the use of high purity cleaning 
solutions (electronic grade solvents and acids) and 
deionized water. Also the metals and chemicals used 
were purified (contact metals 99.99% silver, 99% tita- 
nium, high purity P0C1 3 as the diffusant source) .The 
gases (nitrogen, oxygen, or forming gas) used during 
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diffusion, annealing, or sintering were dried and were 
of good commercial purity. 

To minimize penetration of any impurities left on the 
silicon surface, the temperature- time schedules were 
reduced as much as possible. The right hand column in 
Table 2 indicates the typical processing temperatures 
involved at each process step. 

The holding fixtures, used during the liquid rinsing, 
and the furnace tubes, were carefully cleaned and were 
used specifically for the separate operations. 

Care was taken in handling (minimum pressure from holding 
fixtures or tweezers) to reduce the chance of mechanical 
damage to the shallow PN junction. 

2.2 Measurement Methods used 

The main measurement methods are described. 

2.2.1 Resistivity Measurements 

Four-probe methods (with suitable correction factor^) 
were used to measure the resistivity of the starting 
ingots and slices and for sheet resistance measure- 
ments on the diffused layers. in addition, measure- 
ment of the zero bias capacitance provided a cross- 
check on the dopant concentration at various stages of 
cell processing. 

2.2.2 Etch Pit Density Measurements 

The Sir tl etch (2.1.2 (iv) above) attacked edge disloca- 
tions preferentially, leaving tetrahedral pits on the 
(111) plane. The density of these etch pits was esti- 
mated by photographing known areas at typical locations 
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across the slice and by counting the number of etch 
pits in that area. The results obtained using the Sirtl 
etch are given in Table 3. with a few exceptions, the 
counts are seen to agree with the specifications given 
for the ingots, other etchants (e.g., Dash etch*) were 
used and gave similar densities to the Sirtl etch. 

2.2.3 Solar cell Measurements 

On samples which had been processed to form solar cells, 

measurement of the i-v characteristics was made by varying 

the load resistance from zero to a high value, while the 

cell was illuminated. From the i-v curve, I , V , and 

sc oc 

the maximum power, P , could be read-off; also the CFF 

mdx 

which equals P__ v /l _ x V . in addition, an indication 

Iug X SC o c 

of the spectral response was possible by measurement of 
I under either the xenon or tungsten lamps which were 
blended in the simulator to provide the AMO spectrum. 

Thus, the xenon I showed the response for short wave- 

5C 

lengths (<C 0.6 |un) , whereas the I under tungsten gave 
the response for longer wavelengths (>0.6 yUm) . Detailed 
spectral response curves were measured using a filter 
wheel monochromator. 

2.2.4 Diffusion Length Measurement 

The carrier diffusion lengths (or lifetimes) were mea- 
sured over the wide range of resistivities (0.01 to 10 
ohm-cm) and the consequent expected wide range of dif- 
fusion lengths ( < 10 to > 200 ^tm) . The factors considered 
in choosing the measurement method are given in Appendix A 
Appendix B describes the chosen method (the photovoltaic 
method) , both from theoretical and practical aspects. 

TWo methods were used : 

* Dash etch; 1 HF, 3 HNOj* 8 CH^COOH for 16 hours . 
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(a) Back Current Method 

Here a collecting barrier (diffused pn junction or 
Schottky barrier) was formed on one side of the 
slice together with an ohmic contact? the other 
side of the slice was sandblasted to ensure high 
surface recombination velocity at that surface, and 
an open structure ohmic contact was applied, allowing 
the test illumination (LED, etc.) to fall directly 
on the sandblasted silicon. This structure is shown 
in Figure 5. 

(b) Front current Method 

In this method, the test illumination fell directly 
on the collecting barrier (through an open struc- 
ture contact) the other side of the slice had an 
ohmic contact applied over the slice. 

2.2.5 Diode Measurements 

For no illumination, the diode i-v characteristics could 
be measured on samples which contained a PN junction. 

The i-v characteristics were measured for either forward 
or reverse bias polarity. Also the same dark diode 
characteristics could be measured, as a function of 
position over the whole cell area, by diode measurements 
made on small mesas, formed as in 2.1.3 (v) . 

2.2.6 Trap Density Measurements 

I?he presence of harmful impurities can lead to traps 
which reduce carrier diffusion length. 

Three methods were tried, to find trap densities. 

(a) Variation of Diffusion Length with Temperature 

Careful measurement of L as a function of tempera- 
ture (either increasing or decreasing from room 
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temperature) , with simplifying assumptions about 
the trap properties, could lead to information con- 
cerning trap densities, and their level in the for- 
bidden energy band of silicon. 

(b) Thermally Stimulated Currents 

References (19-21) describe this measurement method. 

A sample with some mesa diodes was cooled to liquid 
nitrogen temperatures. The diodes were illuminated 
with penetrating radiation, thus exciting excess 
carriers into traps, with a constant reverse voltage 
applied to the diodes, the sample was warmed up to 
room temperature. Whenever the temperature provided 
enough thermal energy to empty a trap level, a burst 
of excess dark reverse current was observed. The 
temperature at which the burst occurs could be used 
to show the trap energy level. The magnitude of 
the current pulse at this temperature showed the 
density of traps at that level. 

(c) Thermally Stimulated Capacitance 

Reference 21 describes this method which is similar 
in principle to method (b) , but which offers some 
advantages. 

Again a mesa diode was used. The capacitance at 
large reverse bias was measured for two initial star- 
ting conditions, either by cooling the diode (under 
zero bias) to low temperature (liquid nitrogen), then 
applying reverse bias and warming the diode, or by 
applying the reverse bias at room temperature before 
cooling, and again measuring the change in capaci- 
tance as the sample returned to room temperature. 
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3.0 


EXPERIMENTAL RESULTS 


3.1 Diffusion Length Results 

3.1.1 Diffusion Length versus processing Step 

Figures 6 through 9 summarize the L-values measured on 
fifteen ingots* after the main process steps, the steps 
being shown on the bottom of the figures. An envelope 
curve is used to indicate the range of measured values. 
Generally, the spread in L-values is reasonably low, 
allowing the plotting of average L-values as given in 
Figure 10. On this figure, the ingot parameters have 
been included, it can be seen that the diffusion length 
did not suffer a marked change at any step as was origi- 
nally expected before this program. There did not appear 
to be any correlation between the direction of L-change 
and the resistivity or the oxygen content. The low 
dislocation density ingots showed some tendency to 
decreased L after diffusion. 

3.1.2 Diffusion Length versus Silicon properties 

Using L-values for samples completely processed as cells 
(i.e., after Step VI in Figure 10), L is plotted versus 
resistivity in Figure 11. on this figure, each ingot 
point is identified by two letters, the first letter 
indicating the extreme of oxygen concentration, the 
second letter showing the extreme of dislocation density. 

L fell off rapidly for resistivity below 0.05 ohm-cm 
(for N a > 10 cm ). Also using published values' ’ 
of carrier concentration versus resistivity, l is plotted 
as a function of the acceptor concentration, N , in 
Figure 12. Reference 16 also gives the measured varia- 
tion of carrier mobility y^) with carrier concentration. 

* The two remaining ingots, 973-8 and 9 are discussed separately 
in section 3.3.1 below. 
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Since the diffusion coefficient for carriers is given 
by the Einstein relation 

Vip 

D = (5.) 

and since 

L 2 = DT = (6) 

(where T is the minority carrier lifetime) , 

use of the ^-values corresponding to various values 

results in a set of 7-values corresponding to the l- 

values measured. in this way a plot of T versus 

can be obtained from Figure 12 and is shown in Figure 13. 

Also shown in Figure 13 is a plot of carrier lifetime 
versus the concentration of gold, a widely-used lifetime 
killer. 

3.1.3 General Comments on L-measurements 

L did not change appreciably with the various process 
steps. In general, L decreased as the background doping 
was increased. For any given carrier concentration, 
the L-values for low oxygen concentration ingots were 
higher. However, these low oxygen ingots were grown 
by some modification of float-zone refining, so that 
it could be that the higher diffusion length resulted 
from improved crystal growth conditions rather than 
from the low oxygen content. 

The plot of T versus the gold concentration, in Figure 13 

was included to see if the rapid decreases in ?~with 

18 —3 

dopant concentration for y 10 cm could be attri- 
buted to the accidental incorporation into the silicon 
of an impurity like gold (or iron) , carried in with the 
boron dopant source, if the dopant source contained a 
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fixed percentage of such an impurity (say ^0.1%), then 

part of the 'T-decrease at the increasing N values could 

.A 

be caused by the increasing concentration of the killer 
impurity. 

The dislocation densities did not appear to have a strong 
effect on L, since for each resistivity range comparable 
L-values were obtained with either high or low disloca- 
tion densities. There has been speculation* that since 
almost all the ingots with high dislocation density 
showed increased L after diffusion, there was the possi- 
bility of dislocation-gettering, i.e., the dislocations 
aided in off-setting some of the impurities responsible 
for decreasing L. 

3.2 Solar cell Results 

Solar cell measurements were made for samples processed 

through step VI. From the i-v curve traced under AMO 

2 

illumination, intensity 140 mW/cm , (and the cells at 
28°C), the main photovoltaic parameters, short circuit 
current (I ) and open circuit voltage (v ) were extracted 
and plotted. 

3.2.1 Cell Current Values 

Figure 14 shows how I varied with N . Figure 15 shows 

SC A 

the plot of I versus diffusion length, L. (For compari- 
son. Figure 2 is replotted as a dashed curve on the same 
scale.) The ingot numbers and their resistivity are given 
alongside the experimental points. in many cases, after: 
cell measurements had been made, the back contact was 
partially removed, and the L-value for the cells was 
measured directly by the back current version of the 
photovoltaic method for measuring L. 

* private communication from H. W. Brandhorst, Jr. , NASA-Lewis 
Research Center. 
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As L decreases, I is also seen to decrease. The 

SC 

experimental slope differs from that of Figure 2, 
because that figure was obtained by reducing L (by 
electron irradiation) while maintaining the original 
resistivity (10 ohm-cm) . For the ingots in the present 
study, L was decreased by the additional acceptor doping, 
and the material quantities connecting the diffusion 
length and 1 (mobilities, diffusion coefficients) 
will be different for different doping levels. 

Spectral response curves are plotted in Figures 16 and 

17 along with typical L-values. The losses in I at 

sc 

higher doping levels (lower l) show reduced collection 
(by diffusion) of carriers generated deeper in the sili- 
con by more penetrating wavelengths. 

3.2.2 Cell Voltage Values 

Figure 18 shows V qc as a function of carrier concentration. 

The highest V oc values occurred around N = 7 x 10"^ cm~^ 

A 

(resistivity 0.1 ohm— cm) and fell off predictably at 

lower concentrations, and for concentrations above 10^^ cm~^ 

were less than the simple theoretical v values plotted 

oc . 

in Figure 4. 

The experimental curves in Figure 4 and Figure 18 do not 
agree exactly; the reason is that Figure 4 was based largely 
on measurements on p+/N cells, and since ^ >^, to give 

equivalent resistivities requires concentrations of donors 
srnaller than the acceptor concentration. The data given 
in Figure 18 represents more systematic information for 
P-ingots, over a wider doping range, than was available 
earlier. 

19 


Below (in section 4.0) some speculation is given on the 
possible causes for the reduced V qc values at high doping 
concentrations. 

3.3 other Results 

3.3.1 Measurements on ingots with High L-values 

For two of the ingots (973-8, 973-9) with resistivities 
approximately 10 ohm-cm, the vendor measured very high 
minority carrier lifetimes (400 ^us, 800 respectively, 
corresponding to L ^ 1150 and 1625 yum) . 

These ingots were treated somewhat differently from the 

other ingots, because it was hoped that they would give 

very high I values, and would also allow checking of 
s c 

the assumptions made in the photovoltaic method selected 
to measure L-values. 

Unfortunately, this was not the case. The solar cell 

measurements on slices cut from these special ingots 

did not show any improvement in I resulting from the 

s c 

high L-values; in fact, the spectral response curves for 
cell from these ingots showed reduced long wavelength 
response. 

The L-values measured on both undiffused and diffused 
specimens of these ingots were high, in the range of those 
quoted by the vendor (who used a photoconductive decay 
method for '("-measurement) . However, the various methods 
used to analyze the photovoltaic L-method were not inter- 
nally consistent. More details of the measurements on 
these ingots is given in Appendix B. 

3.3.2 Unilluminated Diode characteristics 

Additional information to explain the differences in 
cell voltage measured for the various ingots was obtained 

20 


by taking typical cells using silicon of 10, 1 , or 0.1 
ohm-era, and dividing the cells into several separate 
mesa diodes (by masking and etching) . The forward biassed 
X-V characteristics with the diode illuminated were 
then plotted for each mesa diode. 

in general it was found that at the higher resistivity 
( ^ 1 ohm— cm) these small area diodes had very uniform 
behavior, but that the spread in characteristics became 
progressively wider as the resistivity decreased.. Figure 19 
shows the measured spread in forward voltage (V J needed 
to produce a forward current equal to 35 mA/cm (the 
operating current density of the diode for an illuminated 
solar cell) . Although the maximum V f increases for the 
lower resistivities, the increasing fraction of the mesas 
with low V - results in seriously decreased V for the 

f Ow 

cell, because in the cell ail those mesa diodes operate 
in parallel, and low voltage for any small area can 
decrease the available V . 

U'- 

Figure 20 shows the envelope of the measured In I f versus 

V curves for three resistivities, 
f 

Some other tests were run wherein the PN junction depth 

was increased to approximately one hundred times that 

[ ' 

used in solar cell PN junctions; in most cases the same 
trend in was observed at the lower resistivities. 

Figures 21 and 22 show the envelopes for two resistivities 
for various junction depths. 

The implications of these measurements are discussed more 
fully in section 4.0. 
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3.3.3 Trap Density Tests 

Three different methods (see 2.2.6) were tried, to give 
an indication of the trap density, or of the approximate 
position of the trap level in the forbidden energy gap. 

None of these methods was satisfactory, for the following 
reasons . 

(a) variation of l with Temperature 

L-measurements were made using the I. method over a 

b 

a range of temperatures of the samples. However, 
serious experimental difficulties were met. First 
at lower temperatures, the back illuminated sample 
was covered with a layer of frost, and effective, 
repeatable optical transmission from the light source 
into the sample was difficult. Second, the absorption 
coefficient of the sources was strongly temperature 
dependent, leading to uncertainties in the "constants" 
to be used in the equations. Attempts were made to 
cancel out these experimental factors by measuring 
L using the ratio of the currents on each side of 
a step in the sample. However, as mentioned in 
Appendix B, the L-vaiue obtained from the step samples 
did not agree with that measured on the separate 
sides of the step, thus this method did not give 
reliable results. 

(b) Thermally Stimulated Currents 

Mesa diodes were formed. However, on cooling down 
they already had large (and unstable) reverse leakage 
currents. Later tests used specially made planar 
passivated diodes formed in typical ingots, but no 
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current bursts were seen on following the procedure 
listed above in 2.2.6 (b) . 

(c) Thermally Stimulated Capacitance 

Experimentally, this method appeared promising, but 
in practice there was hysteresis which prevented 
obtaining repeatable capacitance values. Also the 
trap densities were probably very low for high resis- 
tivity silicon (long diffusion length) . For the 
lower resistivity ingots where the trap density could 
be expected to be higher, the large impurity content 
led to high background capacitance, thereby resulting 
in. the need to measure small percentage changes in* 

capacitance. The tests run did not lead to more 

14 

than an order of magnitude indication (10 traps 

3 

per cm for 1 or 10 ohm-cm silicon) . 

3.3.4 Getter ing Tests 

Because it was realized that the lifetime was determined 
by ingot growth conditions, experimental attempts were 
made to reduce the number of effective traps by subjecting 
typical slices to various "gettering" treatments, in 
these treatments, a liquid phase was maintained at the 
silicon surface for a period long enough for fast-diffusing 
impurities such as iron, gold, or copper to reach the 
surface from the bulk of the slice. The liquid phase 
can act as a sink for these killer-impurities, and on 
cooling the sample, and removing the surface layers, it 
is possible to obtain a slice with reduced trap centers. 
Experimentally, gettering tests were run on both 0.01 and 
0 . 1 ohm-cm ingots, using either aluminum or phosphoro- 
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silicate glass as the liquid phase. The results were 
negative in that the measured diffusion lengths after 
gettering were the same as the L-values measured before 
gettering. 


4.0 


DISCUSSION OF RESULTS 


This study began with the assumption that the quality 
of silicon ingots (as shown by minority carrier life- 
time in particular) was altered as a result of solar 
cell processing. The results showed that this was not 
the case, for the processing methods and materials used 
here. There were relatively small changes in the measured 
values of diffusion length and dislocation density. The 
results can be summarized as follows: 

4.1 Diffusion Length Measurements 

The photovoltaic methods analyzed and tested in this 

work proved effective over the wide resistivity range 

(0.01 to 10 ohm-cm) and over the large range of diffusion 

length values (2 to 1600 /Am) . (The corresponding life- 

-2 

time range is even larger, from 10 to 700 jj.s . ) For 
starting ingots, measurements could be made for resisti- 
vities > 0.05 ohm-cm. The L-values obtained could be 
related directly to the current collection properties 
of solar cells. The measurements were repeatable, and 
the accuracy was better than that obtained in most other 
lifetime or diffusion length measurements. No drastic 
change in the basic methods was needed over the whole 
range of measurements. 

4.2 Diffusion Length Values 

Figures 6 through 9 show that L for a given ingot did 
not vary largely for the cell processing used here. 

This implies that the processing methods did not add 
a significant number of harmful impurities. On the 
other hand, l;;rge L-variations were obtained from ingot 
to ingot. The variation of L with other parameters 
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(Figures 11 and 12) showed that L decreased drastically 

17 -3 

for acceptor concentrations above 5 x 10 cm , and 
that in general ingots grown without use of a crucible 
had higher L-values. 

The cause for decreased L at low resistivities depends 

on starting ingot properties and appears to be related 

to the intentional boron doping. The decrease in L up 
18 -3 

to 10 cm resembles the decreased hole mobility to 

this concentration. However, the large decrease in L 

18 -3 

for concentrations above 10 cm is greater than that 
expected from the mobility decrease caused by impurity 
scattering. 

The curves given in Figure 13 suggested that killer- 
impurities like gold are possibly being added with the 
boron. The negative getter ing tests described in section 
3.3.4 suggest that if such impurities are the prime cause 
for L-decrease, their distribution was not altered much 
by the gettering treatments used. 

The L-values correlated well with the cell I values, 

sc 

and also with the cell spectral response. 

Low Resistivity Silicon Results 

The present work added more information on the variation of 
solar cell properties over a wider range of resistivities. 

(a) Several well-controlled ingots at both 0.1 and 0.01 ohm- 
cm were evaluated. These ingots had a wide variety of 
concentrations of oxygen and dislocation density. 

(b) The diffusion length of these low resistivity ingots 
was measured, and the values obtained could explain the 

I__ fall-off observed in cells made from these ingots. 

.sc ■ ■ 
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(c) Some reasons for the cell v behavior (did not 

increase, as earlier theories predicted, with fall- 
off increasing at lower resistivities) were found 
in detailed study of the large area PN junctions via 
measurements of the dark diode behavior of small 
mesas . 

.3.1 Speculation on V Results in Low Resistivity Silicon 

in section 3.3.2, it was observed that as the bulk resis- 
tivity was decreased, there was an increased fraction of 

. i 

mesas with low V f , even though the maximum v f values were 
higher at lower resistivities. Thus a possible physical 
mechanism for decreased V qc for cells was the increasing 
chance of small areas of the junction contributing excess 
forward current. 

Figure 23 (curve with squares) shows the best V values 

oc 

obtained versus the acceptor concentration. Also shown 

-i ' . 

are two curves showing V qc values obtained by considering 
two different theoretical models. Table 4 lists the 
measured values used in these theoretical estimates and 
gives the derived values which were used to plot the 
curves. 

The two theoretical cases are: 

( 22 ) 

(a) Diode Diffusion Theory 

Saturation current density (j qd ) is given by: 



where n^ is intrinsic carrier density, 

N is acceptor concentration, 

■ *» ■ 
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D id diffusion coefficient for electrons, and 
n 

T is lifetime for electrons in bulk silicon. 

The best measured values of 7* n are shown in Table 4 f 

and these values are combined with the D and N 

n A 

values shown to derive J QD . 

in turn ( v oc ) D i s derived from: 

(V ) = v In 

oc D T J 0D 


( 8 ) 


where v. 


kT 


and 


J = short circuit density taken as 35 mA/cm 
s c 

These V values are shown by circles in Figure 23. 
oc 

(b) Space-Charge Recombination Theory 

Saturation current density is given by: 

.t = 3£UE 


OS 


2 T 


(9) 


where n. =. intrinsic carrier density, 

i 

W = space charge layer width, and 
f = effective carrier lifetime in the space 
charge layer. 

Reference to published results showed that a reasonable 

1 

value for L was £? -rrrr T . Using this estimate, 
e 150 n 

together with measured values of W (obtained from 
cajpacitance measurements), gave the values of J qs in 
Table 4. 


Again, 


(V ) - 2v In T 

' oc's T J 


sc 


( 10 ) 


os 


These values are shown by crosses in Figure 23. 
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ini general, V qc appears to be following the diffusion 
diode theory for low acceptor concentrations. The devia- 
tions from this classical theory are more pronounced 
for N £ 10 cm . Analytically the experimental trend 

suggests an increasing effect of the space charge recom- 

1 7 3 

bination centers for N 2- 10 A cm , until finally for 

! _3 ** 

N. >10 cm , the v values are dominated by the 
A oc 

recombination centers. it is not obvious why the V 

oc 

values can be approximated by choosing T in the space 
charge layer much less than the T value measured in the 
bulk. The J q values derived here agree reasonably well 
with extrapolations at high ( j ) or low (J og ) current 
densities for the In l f versus v f plot. 

As discussed in section 5.0 below, much work remains to 
complete the explanation of the V QC values obtained. 

4.3.2 Other Possible Mechanisms 

The speculation in the previous section assumed two limiting 
cases, namely the diode diffusion theory, and space charge 
recombination theory. However, the precise responsible 
mechanism has not been identified. 

(24) 

Other work funded by NASA-Lewis Research Center has 

led to the conclusion that PN junctions formed in low 
resistivity silicon, have excess leakage caused by chan- 
nelling effects. 

Calculations show that usually much higher doping levels 
than those considered here are necessary for the onset 
of junction tunnelling. However, the measured depletion 

O 

layers were thin (down to 500 A) , and with the increased 
density of recombination centers in the space charge 
layer, tunnelling may be possible. 
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It is possible that the failure of schottky barriers found 
at resistivities below 0.05 ohm-cm is also caused by a 
similar mechanism to that which decreases V qc; . 

4.4 Relevance of Program to objectives 

The objective was to study the effects of processing on 
silicon properties, especially on carrier lifetime. 

The program revealed that conventional solar cell fabri- 
cation technology does not have as marked an effect as 
previously thought on the silicon properties. These 
properties are fixed more by the conditions occurring 
during ingot growth. 

For the ingots with high doping concentrations, a fall- 

off in diffusion length was measured, with corresponding 

decrease in cell current. This current decrease was 

accompanied by a decrease in open circuit voltage at 

18 -3 

the highest doping levels used ( > 10 cm ). 

Thus the original program objective was met? however, the 
experimental work revealed new areas of work required 
to elucidate the behavior of cells made from highly doped 
silicon. 

4.5 implications for Future improvements in Cell Output 

in order to increase cell output further by use of lower 
resistivity silicon, two physical processes must be under 
stood, and if possible, offset. 

(a) The cause of decrease in diffusion length must be 
identified. 

(b) The causes of excess diode currents must be isolated, 
and if possible remedial processing applied to allow 
higher voltages. 
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CONCLUSIONS 


The following conclusions can be drawn: 

(a) Present-day solar cell processes did not have a marked 
effect on the properties of silicon, for ingots covering 
a wide resistivity range (0.01 to 10 ohm- cm) , and 
including extreme values of oxygen content and disloca- 
tion density. 

(b) The properties of solar cells made from those ingots 

indicated that the I behavior could be correlated 

with the measured values of diffusion length, and 

that the V behavior could not be easily explained, 
oc 

(c) The physical cause for reduced L at high doping levels 
was not clearly identified. Possibly the cause is 

a combination of effects caused directly by the higher 
concentrations of intentional dopant (boron in this 
case) , and is also determined by increasing concen- 
trations of "killer” impurities, accidentally incor- 
porated in the growing crystal. 

(d) The physical cause for slower increase of V qc could 
lie in detailed examination of the causes of excess 
diode currents, with possible contributions from the 
same higher concentrations of intentional and uninten- 
tional impurities mentioned in (c) . 

(e) For all but the two ingots of very long diffusion 
length, the photovoltaic method of L-measurement 
described in the appendices proved more effective 
than many other methods used previously. 


6.0 


RECOMMENDAT I ONS 


The following recommendations are made: 

(a) To concentrate on ingots with doping concentrations 

17 -3 

greater than 10 cm , to identify the causes of 
reduced diffusion length, and increased diode current. 

(b) To extend the trap density measurements to see if 
they can be used to aid understanding of the reasons 
for low L and higher I -values . 

(c) To test ingots made with very highly purified dopant 
sources, e.g., use multiple zone-refining of the 
dopant ingot before it is used to grow the single 
crystal. This should be possible for the present 
case, since boron does not redistribute strongly on 
zone refining, whereas possibly the candidate metal 
impurities (like' iron, gold, copper, etc.) will segre- 
gate and can be removed. Possibly also, some tests 

of gettering of the dopant source before its use in 
crystal growth could achieve the same objective. 

(d) To continue work on measuring diffusion length, to 
ensure good correlation of this parameter with mea- 
sured cell values. At present the photovoltaic ^method 
has good promise, as does the surface photovoltage 
method (see Appendix B) . 
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APPENDIX A 


DISCUSSION OF METHODS CONSIDERED FOR MEASURING 
DIFFUSION LENGTH 

The measurement of diffusion length (L) was important 

in this study. Several factors were considered in choosing 

the method for measuring diffusion lengths. 

(a) If possible, the method should be electrodeless to 
minimize alteration of the diffusion length by intro- 
duction of metallic impurities during the metal elec- 
trode formation. if this was not the case, the appli- 
cation of suitable contacts had to be performed at 
low temperatures to avoid the possibility that the 
sample preparation would affect the diffusion length. 
Previous work showed that 600°C was the threshhold 
for thermal stresses, and that temperatures below 

400 °C would be preferable to minimize introduction 
of harmful impurities like iron, gold, etc. 

(b) The method should be applicable to both undiffused 
(starting silicon) and diffused samples. 

(c) The method chosen had to be applicable over the whole 
resistivity range used here (0.01 to 10 ohm-cm) and 
over the expected diffusion length range ( < 10 to 
>250 /Am). It was desirable that a single method 

be used over the whole range to avoid the need to 
correlate methods which measured slightly different 
parameters . 

(d) The method used should provide L-values suited for 
correlation with solar cell i-v measurements. 
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A.l 


Possible Diffusion Length Measurement Methods 
To measure diffusion length (or lifetime) of minority 
carriers, two broad classes of methods are available. 

These classes can be described as transient or steady 
state methods. 

(5—8) 

A. 1.1 Transient Methods v 

in these methods the conductivity of the sample is modu- 
lated (e.g., by injection from a short pulse of light^ 5 ' 6 ^ 
or from a pulsed microwave oscillator ^ 7/ ) , and the rate 
of decay of the increased conductivity can be related to 
the carrier lifetime. in tests made early in this study, 
several difficulties were found for transient methods in 
meeting the factors given above. These difficulties were 
first, that when the carrier lifetime was short, the 
recovery pulse was very short and difficult to measure 
accurately; second, when the silicon resistivity was low 
(<1 ohm-cm generally) it was difficult to raise the 
energy of the input pulse sufficiently to produce a signi- 
ficant increase in sample conductivity without at the same 
time altering the assumptions inherent in these transient 
methods. in addition, at these high light or microwave 
intensities some of the boundary conditions became criti- 
cal and the analytical solutions were complicated. 

It was felt that the problems of the transient methods 
could be overcome by the steady state methods which will 
now be discussed in more detail. 

A.l. 2 Steady-State Methods ^ 

in these methods, excess carriers are injected continuously 
(using light, microwaves,^ electrons, or current-injecting 
barriers) and the steady state conditions are used to 
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indicate the diffusion length of the injected carriers. 

For these steady state methods, fast transients do not 
have to be measured, and in addition, it is easier to 
measure steady currents more accurately over a very wide 
range, thus allowing measurement of a wide range of L~ 
values. Three steady state methods were used in this 
work , the first two being described in less detail than 
the third, which was selected for this study. 

(12 13) 

(a) Surface photovo3.tage Method ' 

For this method, a transparent, conducting contact is 
capacitatively coupled to the silicon and the photovoltages 
induced by light sources of different wavelengths are 
measured. if the equivalent light flux densities to 
produce the same surface photovoltage are plotted versus 
the reciprocal of the absorption coefficients for the 
illuminating wavelengths, the intercept of the line obtained 
on the zero flux density axis gives the diffusion length. 
Some workers have reported good success with this method 
even down to 0.1 ohm-cm. An ohmic contact is required 
for the back surface of the samples, but a low temperature 
contact (pressed indium) has sufficed. initial tests in 
this work did not give repeatable results. However, this 
method has been successfully used both at the National 
Bureau of standards, and also at NASA-Lewis Research 
Center. 

(b) Transistor Gain Method 

Whenever a diffused layer exists on the sample, it is pos- 
sible to form a transistor structure by masking the diffused 
layer on both faces of a silicon slice and etching to 


form an emitter-base-collector transistor structure. 

The current gain for such a wide base transistor can be 
measured and allows extraction of the diffusion length 
in the base. This method was not sensitive where the 
diffusion lengths were much smaller than the base width 
which was generally the case for lower resistivity mate- 
rial. Also for low resistivities there were large para- 
sitic circulating currents in the transistor baseband 
these currents reduced the sensitivity of this method. 

(c) Photovoltaic Methods^ 9 * 

For these methods carriers are provided by a light source 
and the resulting current and voltage at a collecting 
barrier are determined by the diffusion length, if suitable 
allowance is made for the boundary conditions. Sample 
structures as shown in Figure 5 can be used. Three main 
components of a sample can be seen--namely , the collecting 
barrier (a potential barrier which can collect minority 
carriers), a contact to this collecting barrier, and an 
ohmic contact to the silicon sample. 

Two modes of illumination can be used as shown in Figure 5(a) 
and 5(b). For 5(a) (the back current, method) the sample 
is illuminated on the surface opposite to the collecting 
barrier, and the injected carriers must move through the 
slice to reach the barrier and generate the measured photo 
current, in order to obtain measurable currents generated 
in material with very short diffusion lengths, very thin 
slices must be used. in 5(b) (the front current, i f 
method) the sample is illuminated on the side with the 
collecting barrier. The carriers generated are collected 
at the barrier. This mode is more suited to silicon with 





short diffusion lengths, since carriers are generated 
nearer the collecting barrier. For both modes the con- 
tact structure must have open areas on the face which 
is illuminated; typical examples are shown in Figure 5(c). 

The use of these methods to measure diffusion length 
results in structures and measurements which closely 
resemble those used in actual solar cells. Therefore, 
correlation with the measured cell properties is direct. 

The theoretical and practical details of the photovoltaic 
methods are given in Appendix B (below) . 

(d) Electronvoltaic Method 

Another method to inject carriers uses moderate energy 
( '•■' 1 mev) electrons at low injection levels. This method 
requires separate measurement of some parameters (the 
number of hole electron pairs generated per incident 
electron) by subsidiary measurements, and has been regarded 
in previous solar cell work as a standard method. How- 
ever, for the present work, a rapid, easily available 
method was needed. As will be seen in the next section, 
electronvoltaic measurements were made* on selected samples', 
which were then used to calibrate the photovoltaic method. 



* At TRW Systems. 
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DETAILED DISCUSSION OF PHOTOVOLTAIC 
DIFFUSION LENGTH MEASUREMENTS 

B.l Theoretical Basis 

The theoretical basis for the photovoltaic method for 
determining diffusion lengths is summarized here, giving 
the assumptions and practical equations used. The geo- 
metric arrangement is shown in Figure 24, 

B.1.1 List of Symbols Used 

o C - absorption coefficient for photon-source 
D = diffusion coefficient for minority carriers 
T = lifetime of minority carriers 
L = diffusion length of minority carriers 
incident photon flux 
n = density of carriers 
S = surface recombination velocity 
g = generation rate of carriers 
q = electronic charge 
A = sample area 
W = sample thickness 

B.l. 2 General Assumptions Made in Deriving Equations 

(a) Monochromatic illumination was used, characterized 
by a flux (fc 0 and an absorption coefficient oC„ 

(b) One dimensiona^l geometry used, i.e. , the sample has 
large lateral dimensions compared to its thickness, 
and the intensity is uniform over the area illumi- 
nated. 
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(c) unity collection efficiency obtained at the collecting 
barriers (PN junction or metal-semiconductor contact), 
i.e., any minority carrier reaching x = W, will be 
collected and contribute to current. 

(d) A shallow PN junction depth such that 3* •= junction 

1 J 
depth « — (i.e., negligible number of carriers photo- 

generated in diffused layer for i f method. 

(e) Steady-state conditions (i.e., all variables indepen- 
dent of time) . 

(f) Homogeneous sample with no internal electric field. 
Therefore, OL, D, and Twill be independent of x, 
and E (electric field) = 0. 

(g) Low-level injection conditions such that <*, D, and T 
are independent of n, and E - 0. 


(h) Back surface (x = 0) is a perfect sink for minority 
carriers (i.e., the back surface has a perfectly 
collecting metal-semiconductor barrier or the sur- 
face preparation is such that S is infinite) . 


B.1.3 


Basic Equations 

The basic equation for the study of spatial and temporal 
variations of minority carriers in semiconductors is the 
continuity equation. The continuity equation using the 
above-mentioned assumptions is: 


n d^n 
g = ? - d— ^ 

4 dx 


tfiis equation is solved subject to the boundary condi- 
tions that at x = 0, n - 0 (i.e., for the -method, the 
illuminated surface is a perfect sink for minority carriers) 


# 
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leading to a derivation of the spatial distribution of 
carriers, once the spatial minority carrier variation 
is determined, the front and back currents (l f , 1^) can 
be determined as follows. 

B.1.4 Back Current (3^ ) 

The general solution for 1^ is still too complex for 
practical use. Thus the experimental conditions are 
adjusted to allow a simplified equation to be used. 

The experimental conditions are such that: 

a — ♦ oo 


Then, 


(ocl) » 1 

(?) 


i = — 

b otw 




(Bl) 


Here I is the current generated by the total incident 
o 

flux (i = qA^ Q ) . 


Equation (Bl) is useful for a determination of L. Once 

the factor has been determined, e.g. by measurements 

oc 

on independently calibrated samples of known diffusion 
length, the diffusion length can be determined from the 


argument of the 


sWl) 


6) 


function. 


An alternate procedure is to use two samples of the same 

ingot (same L assumed) of different thickness W^, Then 

take the I. -ratio, thus yielding 
b 


bl 

I b2 


W 


Sinh 


© 


sinh 


exp 


<W 


(where — >> 1). 

■ . Il 


42 


The value of L can then be determined by an iterative 
procedure. For the method to be practicable it is 
necessary that the sample thickness be comparable to 
ot larger than the diffusion length; ideally W should 
be greater than two or three diffusion lengths, on the 
other hand if the diffusion length is small, the — factor 
may be so large that the back current becomes too small 
to measure accurately. 

B . 1 . 5 Front Current (i^ ) 

For the front current, it is assumed that the back sur- 
face is still a perfect sink for minority carriers (i.e., 
S B , and the front surface pn junction is a perfect 

collecting barrier, if exp is chosen to be »1, then 
it can be shown that 

T _ _ oc l 

f o 1 + ocl 


For this expression to be useful to determine the dif- 
fusion length it is desirable that the factor XL not be 
too large, for if «L » 1, then i f -i i q and i f is no longer 
a sensitive function of L. using a calibrated sample of 
known diffusion length, the diffusion length of the unknown 
sample can be determined by taking the ratio of the for- 
ward currents, 

= ~2 ’ 1 (B4) 

and solving (B4) for the unknown diffusion length. This 

solution is considerably simplified if <£L, and err, are 

1 2 

both «1. Then 



(B5) 
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One possible problem in using penetrating (low ot) illumi- 
nation so as to satisfy this condition is that of reflec- 
tion from the back surface and bottom contact giving rise 
to multiple passes through the sample. 

B.2 Practical Details of the photovoltaic Method 

This section describes the experimental conditions used 
to apply the theoretical equations above to use the photo 
voltaic (PV) method to measure diffusion lengths. 

B . 2.1 Sample preparation 

Section 2.2 described the sample preparation scheme used 
for all ingots. The thick slice (v 0.10 in.) allowed 
estimates of the diffusion length for the as-received 
ingot. The choice of slice thickness was made to con- 
form to the experimental conditions described in section 
B.l. 

The slice surface conditions were chosen, either to pro- 
vide an effective collecting barrier on one surface or 
to fix the surface recombination velocity on the other 
surface in the range required by the chosen solutions 
to the theoretical equations. ' 

in some cases, the PV method could be used on solar cells 
for the I, method, much of the back contact of the cell 

D 

was removed, and the open silicon surface lapped or sand- 
blasted to increase surface recombination velocity. For 
the if method, the cells could be illuminated at the cell 
front surface, 

B . 2 . 2 Collecting Barriers 

As described above, the diffusion length measurements 
depend on the presence of a barrier which can collect 
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photo-injected minority carrier* effectively. Whenever 
a slice had a PN junction formed by diffusion, this PN 
junction could be used as the collecting barrier. For 
undiffused slices, a series of tests showed that metal- 
silicon Schottky barriers could be used as collecting 
barriers over most of the resistivity range used in this 
work. These tests involved forming a diffused pn junction, 
etching off half of the diffused layer and then applying 
a Schottky barrier to the etched portion, open structure 
contacts were applied to the face opposite the two barriers, 
and the resulting from a small area light spot was 
compared for the diffused barrier and the metal-silicon 
barrier. First moving the spot across the face opposite 
the PN junction gave constant i values showing that the 
diffusion length was constant across the slice. The col- 
lection efficiency of the pn junction was assumed to be 
unity. Generally the Schottky barriers had collection 
efficiencies greater than 0.80, and usually were above 
0.90 (i.e., the value for the Schottky barrier area 
was £*0.9 the for the diffused layer area). The repeati- 
bility of the Schottky barrier formation was shown by 
etching off the barrier, replacing it with a fresh bar- 
rier, and measuring l b ? the collected currents for the 
two barriers showed close agreement. 

It was found that effective Schottky barriers could be 
formed on all 10, 1, or 0.1 ohm-cm samples. For the 
0.01 ohm-cm samples, the Schottky barriers were severely 
shunted, and this prevented repeatable measurements. 

A literature search' ' ' revealed that theory predicted 

and experiments confirmed that Schottky barriers formed 


45 



1 




..-4. . . 


mI 


x. 


on silicon with resistivity less than 0.05 ohm-cm, were 
shunted because of the apparent onset of carrier tunnel- 
ling. other attempts (using moderate temperature lithium 
diffusion or liquid metal barriers) were also unable to 
produce effective collection barriers. For these reasons, 
PV L-measurements could not be made for the 0.01 ohm-cm 
samples before they were diffused. 

B.2.3 Light Sources 

The equations given above in section B.l showed that 
various solutions were possible depending on the absorp- 
tion coefficient, ^nd this parameter in turn was set by 
the choice of the light source. Several different sources 
were used, as follows. 

(a) Light Emitting Diodes (LED) 

LED's had several experimental advantages. They emit 
over a narrow wavelength band, and this band can be adjus- 
ted somewhat by varying the impurity doping in the LED. 

Gallium arsenide LED's were used, emitting around 0.9^/tun. 
Higher intensities could be provided by using an array of 
LED ’ s . 

The LED's used were physically small (the enclosing can 
was ^ 3/16" maximum) , and the emitted light was reason- 
ably well collimated. Thus, when the LED was positioned 
close to the silicon surface, small areas could be illumi- 
nated. The LED wavelengths (^0.9 fjia) were absorbed in 
approximately 2 yum depth in the silicon. The minority 
carriers generated in this 2 ^jn thick sheet could be 
regarded as a thin layer of carriers available for move- 
ment towards the collecting barrier. Despite the fact 
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that most of these carriers would be collected at the 
illuminated surface, (because the surface recombination 
velocity was maintained high intentionally) the resulting 
h> was still measurable except for samples with very 
short diffusion lengths. 

A convenient check on the pv experimental conditions 
was afforded by using two LED's of slightly different 
wavelengths (i = 0.91 /um and 0.94 ^un) mounted on the 
same header with a single glass window of high refrac- 
tive index. Comparison of 1^ values on a sample when 
each LED was switched in separately, gave a rapid check 
if the correct combination of sample conditions (collec- 
ting barrier surface recombination, slice thickness, and 
diffusion length) was being used. 

(b) Sources with low Absorption 

For samples with low diffusion length, the front illumi- 
nation method was used. in this case a source emitting 
long wavelengths (with low absorption) was needed. Tung- 
sten lamps (which are rich in longer wavelengths) were 
combined either with a silicon filter (> V thick), which 
prevented all shorter wavelengths from reaching the samples, 
or with a narrow band pass interference filter with trans- 
mission centered around 1.06yiun (accompanied by a 10 mil 
silicon slice to ensure no short wavelength leakage) . 

For most of the L-measurements used in section 3.0, the 
0.91 pum. LED was used for the method and the tungsten 
lamp with the 1.06 p,m filter for the i f method. 

B.2.4 Measurement Circuit 

To measure the true i gc values even if the sample had 
high shunt leakage, the circuit shown in Figure 25 was 
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developed. A current to voltage converter was used. 

The effective resistance across the sample was 0.01 ohm. 

There was a low voltage drop across the sample equal to 

5 

the measured voltage divided by 10 , the voltage gain of 
the operational amplifier. 


B.2.5 Calibration procedures 

As described above, the boundary conditions on the samples 
made were adjusted to give equations (Bl) and (B'l) for the 
1^ and if respectively. 

These equations were* 


I 

o 

(9 

(Bl) 

” 0CW 

sinh (r) 



(XL 

(B3) 

— X 

o 

(1 + (XL) 


Several different caliJoration procedures were used. 


(a) Use of Samples Calibrated by Electronvoltaic 
Measurements 

Several solar cells covering a range of diffusion lengths 
had their diffusion length measured* by the electronvoltaic 
method (EV method) . 


These calibrated cells were then processed to provide 
samples suited for photovoltaic measurements of diffusion 
length using either front or back, illumination. 

For each series of L-measurements these calibrated samples 

covering the expected L-range were used to check the 
I 

"oc (° r I 0 > value for the particular light source, spacing, 
etc. 

* Measurements made by j. r. Carter, jr., of TRW Systems. 
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(i) Method 

Table 5 lists the EV measurements of diffusion length 
for five cells which had an open back contact to allow 
back current measurements. The cell thicknesses (W) and 

values of the function are given, using equa- 

tion (Bl) I q /o£ was computed, and the resulting values 
are also tabulated, if conditions were maintained con- 
stant, this value of I /at could be used for subsequent 

o 

measurements using the back current, in practice, the 
LED-sample spacing and other conditions varied slightly 
for each series of measurements, but by cross-checking 
with these calibrated EV samples, the applicable I Q /oc 
for each series was derived and used. 

(ii) I £ Method 

Table 6 lists EV measurements of diffusion length for 
five other cells with lower diffusion lengths, on these 
samples, measurements of front current were made, using 
a longer wavelength (lower absorption coefficient) source, 
namely, a tungsten lamp and a 1.06 ^m filter. The i 
values are shown in Table 6, and are plotted versus L in 
Figur i 26. As for the i^-method, for each series of i £ 
measurements, some of these calibrated samples were used 
to derive the constant i in equation (B3) above. 

The absorption coefficient values (ot) in the above equa- 
tions were measured experimentally, as shown in Table 7. 
However, the sources were not strictly monochromatic, 
so that some error was inherent in using "effective 06 - 

values". Where possible, constants such as I /ot were 

o 

derived from the calibrated samples and used in the 
measurements on the various ingots. 
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(b) Use of a Standardized Photocell 

Using a solar cell whose absolute spectral response had 

been measured by the National Bureau of Standards, it 

was possible to estimate the effective photon output of 

the light source around its emitted wavelength, in turn 

the photon current entering the silicon sample could be 

calculated by applying applicable transmission coefficients. 

Figure 27 shows the absolute spectral response of a cell 

used to measure the photon output of two light sources . 

Thbje sources were LED's with emitted wavelengths around 

0.91 ^im and 0.94 these wavelengths being indicated > 

by -the arrows. Table 8 shows the derivation of the photon 

outputs of the two LED's. This method gave I values in 

o 

fair agreement with I q 's calculated from the measure- 
ments, using the absorption coefficient estimated for the 
LED in use. As before, unless the spacings, LED's, etc., 
were maintained constant, it was necessary to recalculate 

the I value for each different series of measurements. 

o 

For comparison, for the 0.91 jam LED, the value of iVoi 

obtained for the EV samples was 5.2 ^-cm as shown in 

Table 5. using the I value derived in Table 8 (890 fx. A) 

° -l 

and the effective ctf-value from Table 7 (254 cm ) gives 

I /oL & 3.5 u A* cm. The errors in L-values measured by the I 
o ' ® 

method caused by this difference in I q /ol were ^5% for 
L ~ 65 / ujn, increasing to ^ 15% for L ~ 110 ^um, and 20% 
for L 130 yuma. 

( c ) use of SPV Method instead of EV Method 
Although preliminary measurements with the SPV method 
did not show good consistency, the promising results 
obtained by NASA-Lewis groups suggest that future work 
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in measuring L-values should re-evaluate the SPV method 
both as the main measurement method, and also for cali- 
bration purposes. 


(d) use of Samples Containing a Step 

These samples have a step formed on the back surface, 

giving two regions of different base widths, W^, W£. Figure 28 

shows such a sample. The back currents measured for the 

two sides of the steps are and i^. 

By (Bl), 


J bl 


Similarly, 


o_ 

oCW, 


&) 


o 

oth 


r W. 


Sinh 


(B6) 


I b2 


<*L 


Sinh( W ll 


The ratio 


•^2 


sinh (^'l 


(B7) 


Usually four iterations were sufficient to derive l in 
this equation. L could "also be obtained from the ratio 
Z b2 

hi 


exp 


by using the fact th?*t the sinh ratio is approximately 

(”r w 2) 


Thus measurement of the current ratio, and 


the step height allows solution for l. in prin- 

ciple, this method offers the most direct way to derive 
L-values which can be used for calibration purposes. 


51 


If the assumption is made that L is constant across the 
sample containing a step (an assumption well justified 
from supporting checks of both 1^ and i f taken across 
the sample) , the solution given is independent of the 
exact photon output and effective absorption coefficient, 
in practice, however, L-values obtained on step samples 
did not agree closely with measurements made using the 
EV samples for calibration. 


An example of such inconsistency is given below. 

A step sample gave the following measurements for the 
two sides of the steps , 


Using 

I 



! bl - 3 ° 

= 528 /m> 

5.2 fji A* cm (from 



X b2 “ 37 A A ' 
W 2 = 422 pja. 


EV samples. Table 5) 



Using the values given, 

sinh 




= 3.27 


W, 


= 2.98 


L = 176 jim. 

Similar calculation for i^, w, gives 

L = 139 fm. 
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The reason for this inconsistency was not found. 

(e) use of ingots with High L-values 

Ingots 973-8 and 9 were used. For L-values £ 1000 yam, 

in equation (b3) , ocl >^> 1 and the measured front current 


(B9) 


in equation (Bl) 
Thus, 

I 


If Iq ’ w 
for small W, ~ is small, and sinh ^"~ 

Ii 


1W 

II 


I 


I q / — 2 . 

b a w 


i.e., I I • oc w 
o b 


(BIO) 


in addition to this limiting condition for reasonably 
thin samples with high L-values, another check was obtained 
by measuring a series of samples of increasing thickness. 
These ingots were sliced to give a series of samples of 
increasing thickness, one group of such samples had a 
Schottky collecting barrier applied, and the other group 
was diffused. Back current L -measurements were made on 
both series of samples; front current L-measurements were 
also made on the diffused group. 


The current measurements made on the samples from ingot 
973-9 are shown in Table 10. These results are now ana- 
lyzed in several different ways. 


(i) Using the 


value obtained from EV calibration 
I 

~ 5 • 2 /<A - cm (from Table 5) 

Table 11 shows the L-values calculated from the t values 

b 

in Table 10. 

(ii) Using equation (BIO) for the thinnest sample. (#1) 

I 

o , „ 

"£* = 1.4 fih 

Using this value for samples 2 through 6 gives the l- 
values shown in Table 12. These values differ from those 
in Table 11 (much higher) . 


(iii) Using measurements for successive samples 

Assuming l is the same for two samples of successive 

thickness, L can be obtained by taking the t ratio for 

b 

the two samples. 


* 

+ 


For example, for samples #4 and #5 in Table 10. 


W 4 

Let T = x 4 . 


w r 


L X 5* 


d4 
■^5 

Assume the 

From (BIO) , 
L 


sinh x 
sinh x 


exp (-2x 4 ) 


fH 


i = g( x 5 ~ x 4 > fl - exp <~ 2 x 5> 
exp (-2x 4 )“] 


I 


1 - exp 


term = 1. 


1074* 


jum 


1698 + yum. 


Iteration 1 using these L-values in equation 
for l gives. 


L = 1140* juun 

refers to undif fused samples, 
refers to diffused samples „ 


2238' 


jum. 


(Bll) 


solving 
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Iteration 2 Using L-values from iteration 1, in (Bll) , 
gives 

L = 1152* jum 2590 + /tm. 

Iteration 3 using L-values from iteration 2, in (Bll), 
gives 


L = 1155* / un 


2860 + jum, 


The values obtained for the undiffused samples are close 
to those shown in Table 12. 

However, the values for the diffused samples differ widely 

The differences between the values Obtained in (iii) and 
(i) are similar to the inconsistencies obtained on the 
step samples ( (d) above) . 

(iv) Graphical Representation 

Equation (Bl) may also be written as (B6) , namely. 


Assume L is constant for all samples. 

I 


Then - k i —p' 

Sinh VLI 


where k. - —— 
1 cCL 


( oc) Assume further that ^ is small. 
L 


Then sinh 


-T 


where = k^L 


I 

o 

06 * 


A plot of In I^ versus In W would be a straight line 
with a slope - -1. Figure 29 shows such a plot for 
the samples of 973-9. For W < 1250 pm, the slope is 


* refers to undiffused samples. 
+ refers to diffused samples. 


close to -1. 
an L-value £ 1000 p.) 


w 

(i.e., since — is small, this implies 


</» 


For the other extreme, assume that — is large. 

fw) L 

fwl 


Then sinh^jH -> £’ 


5.2.6 


Thus (b 6) becomes t = —2 

b 


o 2 
L 


k 3 21 

= “777V where k » 

3 otL 


S'JS 


A plot of In I versus W would give a straight line 
1 D . 

of slope Figure 30 shows the plot. The values 

for L derived from the slope of this plot approaches 

those obtained in (ii) above. This approximation 

has the disadvantage that the i, measurements which 

b 

are most critical for determining L are those for 
the thicker samples where the i^-values are very 
low. 

The fact that back currents were measured (even very low 
values) as a result of minority carriers moving through 

large thicknesses (up to 350 mils) of silicon, supports 
the long lifetime measurements made by the vendor. 

Other Comments on PV Method 

Several other tests were made on the experimental assump~ 
tions underlying the PV method. 

First, transmission comparison made on samples with or 
without the diffused n+ layer, using the low 0C light 
sources used for I ■ -measurements showed that no measurable 
absorption was occurring in the N+ layer. This agrees 
with numerical estimates. Second, transmission tests on 
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silicon slices of varying resistivity showed that the 
effective of-value increased slightly as the resistivity 
of the silicon decreased. These slight ^-variations 
would require adjustments in the calculations if greater 
accuracy in PV L-measurements was required. 

B.3 Summary of the Photovoltaic Method 

The PV method provided a convenient way to characterize 

\ 

a wide variety of ingots. The method showed good repeati- 
bility for samples, and good consistency for several 
samples from the same ingot. Table 9 shows the good 
consistency obtained. This table also shows that there 
was some difference ( ~ 25%) in the L-values obtained by 
the i and i f methods. This is thought to be due to 
slight differences in the analytical assumptions made. 

For higher L-values (> 90 ^m) the j— values of L are 
thought to be more accurate because the assumptions were 
better. For L-values below 90 jx.tr, the i f -values are 
probably more accurate. More work is needed to explain 
the inconsistencies in step samples (d) and long L-ingots (e) 
With these stated limitations, the general level of accu- 
racy compared well with earlier diffusion length or life- 
time measurements, particularly when the large range of 
resistivity and diffusion lengths is considered. 
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Dow Corning 
Dow Corning 
Dow Corning 


Cow 


* All ingots were p-type, boron-doped 
to growth axis. 


TABLE 1 



0.1 > 10 17 > 10 4 

0.1 >10 17 0 

0.1 <2 x 10 15 >3 x 10 4 

0.1 <2 x 10 15 >2 x 10 4 Extra pure 

<10 16 -♦ 0 Additional disloca- 

tion-free pass 

7 >10 17 -'lO 4 

8.5 >10 17 0 

10.5 <2 x 10 15 >2 x 10 4 

17 

9 <10 -%0 Additional disloca 

tion-free pass 

0.01 > 10 17 > 10 4 

0.01 > 10 17 < 10 3 

0.01 < 10 16 — >0 

1.0 > 10 17 > 10 4 

1.0 ?io 17 <io 3 

1.0 < 10 16 -^0 

1.0 <10 16 (?) -#f» 

0.05 ^ <10 16 (?) -*0 

grown in (111) direction and sliced perpendicular 




TABLE 2 


cell processing Sequence 


Step iden- 
tification 
Number 

Operation 

Typical 

Temperature 

°C 

I 

ingot sliced 



Mount slices on plate 
Slices lapped (one side) 



Slices polished (one side) 

-25 


Remove from plate 

(and localized 


heating in 


Clean 

grinding) 

II 

Light chemical etch 



phosphorous diffusion 

— 900 

III 

Anneal 

-860 

IV 

protect polished surface 

Sandblast and/or light etch to remove 
N+ layer on back 

Remove diffusion glass 



Clean for evaporation 

-50 


Evaporate metal contacts (front and bac 

:k) 

V 

Evaporate antireflective coating 

-25 

VI 

Sinter 

— 600 


Edge etch* 

a/ 25 

• 


* Zf processing ended at stages IV or V , this etch was performed. 












TABLE 4 

Saturati on Current Densities (and v \ Calculations 


amp/cm' 


(v ) ' 
oc D 


amp/cm" 


(V )J 
oc'S 


1.1x10" 4.05x10" 


8x10" 


7.3 x 10" 


1.27x10' 


7 x 10 


9 x 10 


2.1 xio" 


5.3 x 10 


615 2.7x10 4.65 x10" 


6 x 10 


6.5 x 10 


1.93 x 10" 


5 x 10 


2.6 x 10* 

6.7 x 10' 


7 x 10*" 3 . 65 x 10*"° 2.23x10' 

5 x 10" 6 5 . 33 x 10“ 10 1 . 08 x 10‘ 


Calculated from resistivity and capacitance measurements. 
Using Irvin's mobility values 16 . 


Calculated from L-measurements. 

J OD = ft 


V oc = 60 lQ gf Q-035 j 

l J OD J 

From capacitance measurements 

Assuming 7~ = —— nr 

1 DL 150 C n 


h. j = 30. i-JL. 

i. V oc = 120 log 












,*,4 . . I . , 4 -w ■ -- . . X i > . . • ,. ...... 4 


TABLE 6 

Electronvoltaic Diffusion Length 
Calibration Samples for i^. Method 


1 

Cell No. 

Thickness 
W, )jLm 

* f + 

(/*&) 

Diffusion 
Length* 
L, /tm 

□ 

oUL 

1+ (OCL) 

I 

o 

(//-A) 

1 - VI - 1 

320 

1,13 

21 

.022 

1.022 

5250 

3 - VI - 3 

311 

242 

73 

.076 

1.076 

3430 

iH 

1 

H 

> 

1 

324 

186 

43 

.045 

1.045 

4319 

5 - VI - 3 

322 

297 

59 

.061 

1.061 

5166 

6 - VI- - 3. 

661 

431 

185 

.192 

1.192 

2676 

Test 15, #6 

305 

281 

102 

.106 

1.106 

2932 


* Measured by electronvoltaic method at TRW Systems 

Forward current pv reading, using tungsten lamp + 1,06 /tun 
filter (+ 10 mils silicon) 

P Calculated using oc= 10.4 cm" 1 . 



63 



_.• i rrr=.*r, ^ r^„ 


TABLE 7 

Measured Absorption coefficient 


Pour Light Sources Used for PV Method 


Light Source 

Effective 4-, 

-1 

cm 

LED, X* 0.91 /tin 

• .... 

254 

LED, X «- 0 . 94 ^m 

154 

Tungsten lamp + 1.06 /*m filter 


(+ 10 mils silicon) 

10.4 

Tungsten lamp + 190 mils silicon 

2 






TABLE 8 


calculation of the Photon Current of an LED 

Using a calibrated photocell* 


. .. — — — 'l 

Quantity 

Light Sources 


LED 1 

LED 2 

Current in LED (mA) 

Wavelength of LED ( A , yLtm) 

Current generated in cell 2A-4B (I gc » 
Cell response at output A (R, fiK/yJH) 
Power output from LED (P) is given by 

40 

0.91 

670 

0.345 

4d 

0.94 

1216 

0.545 

I sc 

Energy per photon at A (E# ©V) 

Number of photons per second output ( <p, sec ) 
p P (watts) 

E E (eV) x 1.6 x 10 19 

1229 

1.36 

15 

5.64 x 10 

2231 

1.316 

16 

1.06 x 10 j 

t 

* 

Current output (1 q ) 



I = q^= 1.6 x 10 (A) 

= 1.6 x 10“ 13 ^ (jaA) 

•! — ■ ■ ■ 

902 jiA 

1694 |iA 

— 1 


* cell 2A-4B was used? its absolute spectral response is shown in Figure 27 
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! TABLE 9 


I and I„ L-measurements on Several Samples 
b ■ 1 

to Show Repeatibility 




Resis- 

W 


■ f 

+ 



ingot 

Sample 

No. 

tivity 

I * 
x b 

L 

*£ 

L 



ohm-cm 

yum 


/4m 

yUA 

/cm 


973-16 

1 

1 

323 

4715 

108 




2 


343 

45.5 

114 

not measured 


undiffused 

3 


327 

45. 

106 





4 


332 

47.5 

111 



973-16 

1 

1 

374 

36. 

116 

840 

67 

Diffused 

2 


386 

36.5 

122 

1099 

125 

3 


389 

37.5 

125 

1057 

114 


4 


379 

40. 

124 

1048 

112 


5 


380 

33. 

115 

1032 

107 


6 


382 

33. 

115 




7 


379 

34.5 

115 

not measured 

.. . 

8 


379 

35. 

116 



; - . , 

973-17 

1 

0.05 

312 

12.5 

64. 



Undiffused 

2 

3 


305 

305 

12.0 

12.5 

61.5 

62. 

not measured 


4 


318 

12.0 

65. 



973-17 

1 

0.05 


2.5 


7 09 

48.6 

Diffused 

2 


368 

5.2 

64. 

778 

57.5 

3 


368 

4.6 

62.5 

742 

52.5 


4 


365 

4.3 

61. 

756 

54.5 


5 


368 

4.2 

61.5 

763 

55.5 


6 


365 

4.5 

61.5 

770 

56.5 


7 


370 

4.2 

61.5 

754 

54.3 


8 


375 

4.0 

62. 

765 

55.7 


9 


368 

4.5 

62. 

762 

55.5 

j 

10 


378 

3.8 

60.5 

726 

50.5 


* Using LED# A' 0.91 /im. 
+ Using LED, A v 0.94 ^m. 


hiMte if. iii diwfl uml - — - ' 
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TABLE 10 

1^ Measurements for Thick Samples from 
Ingot 973-9 (Long Diffusion Length) 


Sample 

undiffused Samples 

Diffused Samples 

No. 

Thickness 

I * 
b 

Thickness 

V 


W, yum 


W, jm 

/4A 

1 

389 

36 

326 

44.6 

2 

653 

20 

660 

23 

3 

1207 

10 

j 1176 

12.4 

4 

1684 

6.2 

1636 

7.2 

5 

2210 

3.8 

2154 

5.3 

6 

8800 

0.1 

8820 

0.09 


Using light source, LED, hr 0.91 jubbl. 


* 









TABLE 11 


Calculated Diffusion Lengths in Thick Samples 
from ingot 973-9 


Sample 

Number 

Thickness 

jusa 

Undiffused Samples 

Diffused Samples 

I b 

- 

L* 

Aim 

U A 

L* 

Ussy 

1 

389 

36 

123 

44.6 

— 

104 

2 

653 

20 

195 

23 

218 

3 

1207 

10 

' 

355 

12.4 

380 

4 

1684 

6*2 

475 

» 

fO 

476 

5 

2210 

3.8 

575 

5.3 

625 

6 

8800 

. ■ ■ . , ■■ ■ 

0.1 

2330 (?) 

0.09 

1300 


Notes used LED with X = 0.91 ^un. 

* Calculated using I o - 5.2 j^.cm, 

cL 
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TABLE 12 

Calculate d Diffusion Lengths in Thick Samples 
from ingot 973-9 


Sample 

Number 

Thickness 

ytcm 

Undiffused Samples 

Diffused 

Samples 


I b 

L* 

r b 

L* 



/UA 

jm 

/*A 

/m 

1 

389 

36 

7700 ( ?) 

44.6 


2 

653 

20 

1010 

23 

— 

3 

1207 

10 

1270 

. 

12.4 

— — 

4 

1684 

6.2 

1220 

7.2 

1560 

5 

2210 

3.8 

1190 

5.3 

1930 

6 

r ■ ■ J, 

8800 

0.1 

1740 

0.09 

1680 


Notes Used LED with A = 0.91^ 

I 

* Calculated using = 1.4 ixA-cm. 
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Open circuit Voltage (mV) 


FIGURE 4 

V vs. Resistivity 

—QC ’ 1 " 

Comparison of Theory and Experiment 

(Ref. 2) 



Resistivity (ohm- cm) 


10 


100 



FIGURE 5 

Sample Structures for py L-measurements 


^ Light 




17(c) Open Area Contact Structures Used 









FIGURE 7 

Diffusion Length after Processing Ste| 
ingots 973- (4 to 7) 
and 973-15 

Ingot No. 



973-6 

973-15 


'3-5 1 


973-4 



FIGURE 8 



i ii hi iv y vi 

Processing Steps 





FIGURE 10 




LEGEND: I - after polishing IV - after contacts 

II - after diffusion V - after coating 

III - after annealing VI - after sintering 


NOTE: ingots 973-8, 9 have L 7 1000 /am at steps I and VI. 



(see legend) 




LEGEND; The first letter - high or low O ^ 

The second letter - high or low dislocation 



it £ 








Short Circuit Current, AMO (mA) 




LEGEND: 


FIGURE 15 

t vs. Diffusion Le ngth 
. “SC ' 1 

0-10 ohm- cm (from Figure 1) 
o - 10 ohm-cm 
+ - 1 ohm-cm 


x - 0.1 ohm-cm 
A - 0.01 ohm-cm 
V - 0.05 ohm-cm 




Spectral Response 





Spectral Response (^A^w-W) 







Open Circuit Voltage (mV) 






Number of Mesas 


FIGURE 19 



























FIGURE 24 

Geometric Arrangement for Theoretical. Analysis 


Photon Flux ( 4> ) for i .-method 
' x o f 



Photon Flux (d> ) for I -method 

1 0 D 













Response (^A/^W) 


FIGURE 27 






FIGURE 28 





for 973-9 (Long L- ingot) 


LEGEND: 


x - undiffused samples 
o - diffused samples 



Back Current 


FIGURE 3Q 

g (Back Current) vs. Slice Thickness 
for 973-9 (Long L-ingot) 

LEGEND: x - undiffused samples 

o - diffused samples 
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Slice Thickness (^un) 




